Strong and weak clustering of inertial particles in turbulent flows 
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We suggested a theory of clustering of inertial particles 
advected by a turbulent velocity field caused by an instability 
of their spatial distribution. The reason of the clustering in- 
stability is a combined effect of the particle inertia and finite 
correlation time of the velocity field. The crucial parame- 
ter for the instability is a size of the particles. The critical 
size is estimated for a strong clustering (with a finite frac- 
tion of particles in clusters) associated with the growth of the 
mean absolute value of the particles number density and for a 
weak clustering associated with the growth of the second and 
higher moments. A nonlinear mechanism for a saturation of 
the clustering instability (particles collisions in the clusters) 
is suggested. Applications of the analyzed effects to the dy- 
namics of aerosols and droplets in the turbulent atmosphere 
are discussed. The critical size of atmospheric aerosols and 
droplets in clustering is of the order of (20 — 30)fim, and a 
lower estimate of the number of particles in a cluster is about 
hundreds. 



I. INTRODUCTION 

Formation and evolution of aerosols and droplets in- 
homogeneities (clusters) are of fundamental significance 
in many areas of environmental sciences, physics of the 
atmosphere and meteorology (e.g. smog and fog forma- 
tion, rain formation), transport and mixing in industrial 
turbulent flows (like spray drying, pulverized-coal-fired 
furnaces, cyclone dust separation, abrasive water-jet cut- 
ting) and in turbulent combustion, see e.g. [EJ and refer- 
ences therein. Analysis of experimental data shows that 
spatial distributions of droplets in clouds are strongly 
inhomogeneous [0 . Small-scale inhomogeneities in parti- 
cles distribution were observed also in laboratory turbu- 
lent flows |,§. 

It is well-known that turbulence results in a relaxation 
of inhomogeneities of concentration due to turbulent dif- 
fusion Q, whereas the opposite process, e.g., a preferen- 
tial concentration (clustering) of droplets and particles in 
turbulent fluid flow still remains poorly understood. 

In this Letter we suggest a theory of clustering of parti- 
cles and droplets in turbulent flows. The clusters of par- 
ticles are formed due to an instability of their spatial dis- 
tribution suggested in ref. Q and caused by a combined 



effect of a particle inertia and a finite velocity correlation 
time. Particles inside turbulent eddies are carried out 
to the boundary regions between them by inertial forces. 
This mechanism of the preferential concentration acts in 
all scales of turbulence, increasing toward small scales. 
An opposite process, a relaxation of clusters is caused by 
a scale-dependent turbulent diffusion. The turbulent dif- 
fusion decreases towards to smaller scales. Therefore, the 
clustering instability dominates in the Kolmogorov inner 
scale r), which separates inertial and viscous scales. Expo- 
nential growth of the number of particles in the clusters 
is saturated by their collisions. 

II. QUALITATIVE ANALYSIS: INSTABILITY 
GROWTH RATES FOR STRONG AND WEAK 
CLUSTERING 

For an estimate of cluster growth rates we use the equa- 
tion for the number density n(t, r) of particles advected 
by a turbulent velocity field u(t,r): 

dn/dt + V • (nv) = DAn, (1) 

where D is the coefficient of molecular (Brownian) diffu- 
sion. Due to inertia of particles their velocity v(t,r) ^ 
u(t,r). Equation ([j]) conserves the total number of par- 
ticles. Equation for 8(i, r) = n(t, r) — n follows from (Q): 

dQ/dt+ (v ■ V)6 = -Qdivv + DAG, (2) 

where n is the uniform mean number density of particles, 
and the mean particles velocity is zero. We neglected the 
term oc n div v describing an effect of an external source 
of fluctuations which is usually weaker than the effect of 
self-excitation of fluctuations. In order to elucidate par- 
ticles clustering we analyze a role of different terms in 
eq. ^ using a reference frame moving with a cluster ve- 
locity V c \(t). The advective term, [(v — V c \) ■ V]0, causes 
turbulent diffusion inside the cluster with the coefficient 
D T ~ £ c \v c i/3, where £ c \ is the characteristic size of a clus- 
ter and v c \ is the turbulent velocity at the scale £ c \ and 
8(i,r) are the fluctuations of the number density in the 
comoving frame. The term oc — 8 div (v — V c \) can result 
in an instability. Indeed, neglecting diffusion we find the 
solution of eq. (|): ®(t, r) ~ 8(0, r) exp[-I(t)}. Here 

I(t)= / 6(r)rfr = Vj„, I n = / b(r)dr, 

JO n=1 Jnr v 
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and b(r) = div[«(r,r) - Voi(r)] and TV = t/r„ > 1. 
We neglected a spatial dependence of 6(r, r) inside the 
cluster and consider b(r) as a random process with a cor- 
relation time t„ which is of the order of turnover time of 
£ c i-scale eddies, t v ~ £ c i/v c \. Now the integrals /„ can be 
considered as independent random values and the sum 
I(t) is estimated by the central limit theorem: I(t) ~ 
\J (T v b 2 ) v t f, where (...}« is averaging over turbulent ve- 
locity ensemble, ^ is a Gaussian random variable with 
zero mean and unit variance. Now, averaging over statis- 
tics of £ by the identity (exp(fc£))c = exp(fc 2 /2) (with 
k = y/(T v b 2 ) v t ), we estimate (|6| 9 ) ~ (|0 O | 9 ) exp^t) 
with the growth rate j q of the g-th moment given by 

7 ? ~5<T„(div«) 2 )„<z 2 -D T q/£ 2 cl . (3) 

Here the turbulent diffusion inside the cluster is crudely 
taken into account by the the term oc D T . Clearly the in- 
stability is caused by a nonzero value of (^(divu) 2 ), i.e., 
by a compressibility of the particle velocity field v(t, r). 

Compressibility of fluid velocity itself u(t, r) (including 
atmospheric turbulence) is often negligible and divit = 0. 
However, due to the effect of particle inertia their velocity 
v(t, r) does not coincide with u(t, r) ||, and a degree of 
compressibility of the field v(t,r), <7„, defined by 

a v = ([divv} 2 ) v /(\Vxv\ 2 ) v , (4) 

may be of the order of 1 For inertial par- 

ticles divv ~ TpAP/p, where t p is the Stokes time, 
r p = rn p / '6ir pva — 2p v a 2 /Qpv, m p , p p and a are the mass, 
density and radius of particles, respectively. The fluid 
flow parameters are: the viscosity v, density p, pressure 
P, Reynolds number Tie = Lu T /v and the dissipative 
scale of turbulence rj = L7& _3 / 4 . Now we can evaluate: 

°~v - {Pv/pfiahY = (a/a,) 4 (5) 

(see Q), where a* is a characteristic radius of parti- 
cles for which a~ v = 1. For water droplets in the at- 
mosphere Pp/p — 10 3 and a* ~ ?7/30. For the typical 
value of rj ~ 1mm this yields a» ~ 30/im. On windy 
days when r\ decreases, the value of a* correspondingly 
becomes smaller. 

Then we estimate (t v [<Hvv] 2 ) v as 2a v /r v . Taking the 
cluster size £ c \ of the order of the inner scale of turbu- 
lence, 77, we have to identify t v with a turnover time 
of eddies in the inner scale 77, t v — ► = i]/v n ~ 
(L/u T ) lie ^l 2 . Thus, the growth rate 7, may be evalu- 
ated as 

lq — 7ci q{q - ?«) , 7d ~ o- v /t v , q a - l/3a v . (6) 

It is clear that moments with q > q CT are unstable. It 
follows from eqs. (||) and (JsJ) that this happens when 
a > a<;,cr where a 9 , cr = ai^/q 1 / 4 is the value of a at 
which q CI = q. The largest value of a 9jCr corresponds to 



the instability of the first moment, (|0|): ai, C r ~ 0.8a*, 
a 2 ,cr ~ 0.84ai, cr , a 3 ^ cr w 0.76ai jCr , a 4cr w 0.71ai )Cr , eic. 

Note that if (|0|) grows in time then almost all parti- 
cles can be accumulated inside the clusters (if we neglect 
a nonlinear saturation of such growth). We define this 
case as a strong clustering. On the other hand, if g cr > 1 
the first moment (|0|) does not grow and the clusters 
contain a small fraction of the total number of particles. 
This does not mean that the instability of higher mo- 
ments is not important. Thus, e.g., probability of binary 
particles collisions is proportional to the square of their 
density (n 2 ) — [n) 2 + (|0| 2 ). Therefore the growth of the 
2nd moment, (|0| 2 ), (which we define as a weak cluster- 
ing) results in that binary collisions occur mainly between 
particles inside the clusters. The latter can be important 
in coagulation of droplets in atmospheric clouds whereby 
the collisions between droplets play a crucial role in a 
rain formation. The growth of the g-th moment, (|0| 9 ), 
results in that g-particles collisions occur mainly between 
particles inside the cluster. The growth of the negative 
moments of particles number density (possibly associated 
with formation of voids and cellular structures) was dis- 
cussed in H (see also |To|Jll| ). 

III. GROWTH OF THE 2ND MOMENT 

In a previous section we estimated the growth rates 
of all moments (|0| 9 ). Here we present the results of a 
rigorous analysis for the evolution of the two-point 2nd 
moment 

$(t,R) = {e(t,r)Q(t,r + R)) . (7) 

In this analysis we used stochastic calculas (e.g., Wiener 
path integral representation of the solution of the Cauchy 
problem for eq. (0), Feynman-Kac formula and Cameron- 
Martin-Girsanov theorem). We showed that a finite cor- 
relation time of a turbulent velocity plays a crucial role 
for the clustering instability. Notably, an equation for the 
second moment $(t, R) of the number density of inertial 
particles comprises spatial derivatives of high orders due 
to a non-local nature of turbulent transport of inertial 
particles in a random velocity field with a finite corre- 
lation time ||. However, we found that at least for two 
models of a random velocity field [(i) the random velocity 
with Gaussian statistics of the integrals J Q v(t' , £) dt' and 

f b(t',£)dt'; and (ii) the Gaussian velocity field with a 
small yet finite correlation time] the equation for Q(t, R) 
is a second-order partial differential equation: 

d$/dt = £$(t,R), (8) 
t = B{R) + 2U{R) ■ V + D a(3 (R)V a V(3 • 

For the first model of the turbulent velocity field, e.g., 
the coefficients in eq. (||) are given by: 
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B(R)*2 (6[0,C(ri|0)]6[T^(r 2 |T)])(iT, (9) 
Jo 

poo 

U(R)n-2 / (v[0,t(ri\0)]b[T,e(r 2 \T)])dT, 
Jo 

T />OG 

D af3 (R)^2 / (watO.fCnlOj^tr^Cralr)])^, 
Jo 



where D a p(R) 

~ T „ T 



2D6 a(i + D ap (R) and £> Q/3 (i?) = 

D a p(0) — D a p{R) is the scale-dependent tensor of tur- 
bulent diffusion. 

For the ^-correlated in time random Gaussian com- 
pressible velocity field the operator C is replaced by Cq 
in the equation for the second moment R), where 

C = B (R) + 2U (R) ■ V + D aP (R)V a y p , (10) 
B (R) = \7 a VpD aP {R), U o , a (R)=V D a0 {R) (11) 



(for details see |p[|l2||). In the (^-correlated in time velocity 
field the second moment <£>(i, R) can only decay in spite 
of the compressibility of the velocity field. The reason 
is that the differential operator C = V Q.V pD a p(R) is 
adjoint to the operator Cq = D Q| a(i2)V a V 'p and their 
eigenvalues are equal. The damping rate for the equation 



d$/dt = £mt,R 



(12) 



has been found in ref. |T^] for a compressible isotropic 
homogeneous turbulence in a dissipative range: 
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6t v (1 + er T )(l + 3cr T ) ' 



(13) 



where a T is the degree of compressibility of the tensor 

T 

D a p{R). For the (^-correlated in time incompressible ve- 
locity field (cr T = 0) eq. ( |l2| ) was derived in ref. Q. 
Thus, for the model of turbulent advection with a delta 
correlated in time velocity field the clustering instability 
of the 2nd moment does not occur. 

A general form of the turbulent diffusion tensor in a 
dissipative range is given by 

D T a/3 (R) = (dR 2 S aP + C 2 R a Rp)/T v , (14) 
Ci = (4 + 2a T )/3(l + a T ), C 2 = (4cr T - 2)/3 (1 + a T ). 

The parameter <t t is defined by analogy with eq. (0) : 
V D T V V a V D T a0 (R) 



> (15) 



where e Q /3 7 is the fully antisymmetric unit tensor. By def- 
initions (Q) and ( |l5| ) a T = a v in the case of 5-correlated in 
time compressible velocity field. Equations (^) show that 
for a finite correlation time identities (^) are violated: 
B(R) + B (R), U(R) ^ U (R). 

Let us study the clustering instability. Equation (^) in 
a nondimensional form reads: 



~dt 



m(r) Lm(r) 



+ (f/ - C 2 )r 2 



2$' 

r 



l/m(r) = (Ci+C 2 )r 2 + 2/Sc 



(16) 



where U = UR and Sc = j//£) is the Schmidt number. 
For inertial particles Sc 3> 1. The nondimensional vari- 
ables in eq. ( |l6| ) are r = R/rj and t — t/T v , B and U are 
measured in the units tZ . In a molecular diffusion re- 
gion of scales whereby r <C Sc -1 / 2 , all terms oc r 2 (with 
Ci , C 2 and U) may be neglected. Then the solution 
of eq. ( ^6|) is given by <&(r) = (1 — ar 2 ) exp(7 2 i), where 

a = Sc(B — 7 2 T,j)/12 an( i B > r )2 T r r Id a turbulent dif- 

i / 2 

fusion region of scales, Sc ' < r C 1, the molecular 
diffusion term oc 1/Sc is negligible. Thus, the solution 
of eq. ( [l6| ) in this region is <&(r) = A\r~ x exp(7 2 i), 
where A = {C 1 - C 2 + 2U ± iC 3 )/2(Ci + C* 2 ), and 
C| = 4(5 - 7ar,)(Ci + C 2 ) - (C x - C 2 + 2C/) 2 . Since 
the total number of particles is conserved in a closed vol- 
ume, J °° r 2 <I>(7-) dr — 0. This implies that C 2 > 0, and 
therefore A is a complex number. In addition, the corre- 
lation function $(r) has global maximum at r — 0. This 
implies that C\ > C 2 — 2U. The latter condition, e.g., 
for very small U yields <tt < 3. For r » 1 the solution 
for <&(r) decays sharply with r. The growth rate 7 2 of 
the second moment of particles number density can be 
obtained by matching the correlation function <f>(r) and 
its first derivative $'(r) at the boundaries of the above 

-I Icy 

regions, i.e., at the points r — Sc ' and r = 1. The 
matching yields C 3 /2(d + C 2 ) « 27r/lnSc. Thus, 
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200^0^-^) (3-<r T ) 5 



r,,(l + 3cr T ) 
3tt 2 (1 + 3(7, 



3(1 



(1 + cr T )ln Sc 



20(a E 



6(1 + CT T ) 

•o"u) 



(17) 



t v (1 + <r B )(l + <r v ) 



where we introduced <r B and a v defined by equations B = 
20ct b /(1 + <j b ) and U = 20ct u /3(1 + a v ). Note that the 
parameters er B w a u ^ a v . For the 5-correlated in time 
random compressible velocity field o~ B — o~ v = a T = a~ v . 
Note that eq. (|l3|) is written for Sc — > oo. Analysis of 
eq. ( p7| ) shows that the critical value of o~ v required for 
the clustering instability is cr cr w (0.1 — 0.2). Notably, 
in the second model of a random velocity field (i.e., the 
Gaussian velocity field with a small yet finite correlation 
time) the clustering instability occurs when o~ v > 0.2. 
Equation (^|) also yields for the threshold of the insta- 
bility of the 2nd moment (at g cr = 2) a similar value 
cr cr ~ 1/6. 



IV. NONLINEAR EFFECTS 

The compressibility of the turbulent velocity field with 
a finite correlation time can cause the exponential growth 



3 



of the moments of particles number density. This small- 
scale instability results in formation of strong inhomo- 
geneities (clusters) in particles spatial distributions. The 
linear analysis does not allow to determine a mechanism 
of saturation of the clustering instability. As can be seen 
from eq. ( |l7| ) molecular diffusion only depletes the growth 
rates of the clustering instability at the linear stage (con- 
trary to the instability discussed in rcf. || ) . The cluster- 
ing instability is saturated by nonlinear effects. 

Consider now a mechanism of the nonlinear satura- 
tion of the clustering instability using on the example of 
atmospheric turbulence with characteristic parameters: 
t] ~ 1mm, ~ (0.1 — 0.01)s. A momentum coupling of 
particles and turbulent fluid may be essential when the 
mass loading is not small: m v n c \ ~ p. For this condition 
the kinetic energy of fluid p(u 2 ) is of the order of the 
particles energy m p n c \(v 2 ), where \u\ ~ |t?|. This yields: 

n c \ ~ a~ 3 (p/3p p ) . (18) 

For water droplets p p /p ~ 10 3 . Thus, for a = a* ~ 
30/im we obtain n c \ ~ 10 cm -3 and the total number 
of particles in the cluster of size ry, N c \ ~ 77 3 n c i ~ 10. 
These values may be considered as a lower estimate for 
coupling effects of particles on the flow. 

An actual mechanism of the nonlinear saturation of 
the clustering instability is associated with particles col- 
lisions causing particles effective pressure which prevents 
further grows of concentration. Particles collisions play 
essential role when during the life-time of a cluster the 
total number of collisions is of the order of number of par- 
ticles in the cluster. The rate of collisions J ~ nci/r-q can 
be estimated as J ~ 47ra 2 n 2 1 |'Urei|- The relative velocity 
v Te \ of colliding particles with different but comparable 
sizes can be estimated as |w rc i| ~ t p \(u ■ V)tt| ~ T p u 2 /n. 
Thus the collisions in clusters may be essential for 

n cl ~a- 3 ( V /a)(p/3p p ), £ s ~ a{iap v /r,p) 1 ' 3 , (19) 

where £ s is a mean separation of particles in the cluster. 
For the above parameters (a = 30/im) n c i ~ 3x 10 5 cm~ 3 , 
4 ~ 5a « 150/mi and N c \ ~ 300. Note that the mean 
number density of droplets in clouds n is about 10 2 — 
10 3 cm~ 3 . Therefore the clustering instability of droplets 
in the clouds increases their concentrations in the clusters 
by the orders of magnitude. 

In all our analysis we have neglected an effect of sed- 
imentation of particles in gravity field which is essential 
for particles of the size a > 100/im. Taking t c \ ~ r\ 
we assumed implicitly that r p < t v . This is valid (for 
the atmospheric conditions) if a < 60^m. Otherwise the 
cluster size can be estimated as l c \ ~ r\(r p jT v f'l' 1 . 

Our estimates support the suggestion that the clus- 
tering instability serves as a preliminary stage for a co- 
agulation of water droplets in clouds leading to a rain 
formation. 



V. SUMMARY 

• We showed that the physical reason for the clustering 
instability in spatial distribution of concentration of par- 
ticles in turbulent flows is a combined effect of the inertia 
of particles leading to a compressibility of the particle ve- 
locity field v(t, r) and finite velocity correlation time. 

• The clustering instability may result in a strong clus- 
tering in which a finite fraction of particles are accumu- 
lated in the clusters and a weak clustering when a finite 
fraction of particle collisions occurs in the clusters. 

• The crucial parameter for the clustering instability 
is a radius of the particles a. The instability criterion is 
a > a CI « a* for which ((divu) 2 ) = (|rotv| 2 ). For the 
droplets in the atmosphere a* ~ 30/im. The growth rate 
of the clustering instability 7 c i ~ r^ 1 (a/a*) 4 , where t v 
is the turnover time in the viscous scales of turbulence. 

• We introduced a new concept of compressibility of 
the turbulent diffusion tensor caused by a finite corre- 
lation time of an incompressible velocity field. For this 
model of velocity field, the field of Lagrangian trajecto- 
ries is not divergence-free. 

• We suggested a mechanism of saturation of the clus- 
tering instability - particle collisions in the clusters. An 
evaluated nonlinear level of the saturation of the droplets 
number density in clouds exceeds by the orders of mag- 
nitude their mean number density. 
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